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1. INTRODUCTION 

The effective leptonic weak mixing angle, 
sin^ ^eff** I is a crucial observable in the indirect 
determination of the Higgs boson mass from pre- 
cision experiments up to the Z boson mass scale. 
Defined through the vertex form factors of the Z 
at its mass shell 



sin^^^r = i(l 



3? 



gy(M|) 
<?a(M|) 



(1) 



it is an UV and IR finite quantity. The current 
experimental precision requires a careful analy- 
sis of the error on the theory side. Here, sev- 
eral contributions have been calculated over the 
years. In particular, two- and three-loop QCD 
|1I2| . two-loop electroweak fermionic 3 4 , and 
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parts of the three-loop electroweak contributions 
|S] are known. Recently, also the Higgs bo- 
son mass dependence of the purely bosonic elec- 
troweak graphs has been obtained 0. The com- 
putation of the W boson mass from fi decay, 
which is part of the final prediction for sin^ ^Iff * 
is more complete, since the two-loop electroweak 
corrections have been fully evaluated |7I8| . 

It is important to note, that the prediction from 
IHI has been used by the LEP Electroweak Work- 
ing Group for their final report 9 . It is, there- 
fore, necessary, to make sure that the theory error 
is not underestimated due to the lacking bosonic 
electroweak corrections. To this end, we have cal- 
culated the appropriate diagrams |1(J| and shown 
that they fall well within our estimates. It is the 
purpose of this work to present the details of our 
computation. 

2. MASS DIFFERENCE EXPANSION 

The most complicated bosonic diagrams 
needed for sin^ ^eff * two-loop vertices with up 
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Figure 1. Ultrasoft graphs in the mass difference 
expansion. Dashed hnes go on-threshold, when 
Mw = Mz, and thin sohd hnes are massless. 



to three different masses, Mw , Mz and Mh- It 
is only natural to exploit the small differ- 
ence between Mw and Mz to reduce the number 
of scales. Furthermore, since the authors of [HI 
have evaluated the Higgs boson mass dependence 
of the result normalized at Mh = 100 GeV, it is 
only necessary to evaluate the graphs at the same 
point at first (of course an independent check of 
the mass dependence is also required). This al- 
lows a further mass difference expansion, with ex- 
pansion parameter = 1 — Mfj IM\ . 

It is clear that if there are thresholds at — 
Mz, then starting from some order in the expan- 
sion in Sw = 1 — My^/M'^, we are going to en- 
counter divergences. A suitable technique to re- 
cover the correct value is the expansion by regions 
method ^2]- The idea is to analyze the momen- 
tum regions, which can contribute to the inte- 
gral, and expand the integrand in each region suit- 
ably performing the integration in dimensional 
regularization. In the case, when only one line can 
go on-threshold, there are just two regions: ul- 
trasoft, (us), where fci_2 ^ s^Mz] and hard, (h), 
where fci^2 ~ Mz- The classification of the size of 
the momenta has to be supplied by a momentum 
routing, which clearly needs to be such that the 
selected line goes on-threshold when sw — 0. 

In the (us) region, one encounters three possi- 
ble topologies. Fig. ^ The expansion transforms 
a Feynman integral into an unusual form, since 
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The last integral in Fig.Q] /s, with the above sub- 
stitutions (actually just Eq. Q) has been evalu- 
ated in The other two turn out to be re- 
ducible by the integration-by-parts technique to 
the last case. In practice this exercise has been 
left to IdSolver ^S], with the result 



and 



1 4e-3 ^ 
s2,M| 1 - 2e^ 
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(4) 
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where 



h = 



r2(l-e)r(4e-3) (6) 



X {M'zr''i-s'wr''- 

Since there is no problem with the evaluation of 
the master integrals, the rest of the calculation 
has been reduced to the application of IdSolver 
to all integrals with numerators and dots. 

3. MASTER INTEGRALS 

The integrals in the hard momentum re- 
gion, (h), have the usual Feynman integrand 
structure, with some propagators on-threshold. 
In practice this means that each one of them is 
given by a series in e with numeric coefficients 
multiplied by a trivial scaling factor. The prob- 
lem is now that we have spurious poles in the 
reduction to master integrals and we need ana- 
lytic results for constant parts of some nontrivial 
integrals in order to check the exact cancellation 
of divergences. Of course, one could think of con- 
structing an e-finite basis in the sense of JJ, , but 
we would in the end have to evaluate integrals 
with collinear singularities (even though just to 
finite parts). Therefore we tried to find a basis 
that would produce a relatively small number of 
spurious poles, but in front of simple integrals. 

Indeed, there are 73 master integrals, 1 needed 
to ©(e^), 6 to 0{e^), and 26 to 0{e). Most of the 
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Figure 2. Some hard momentum region graphs 
needed to 0{e). Thick hnes are massive and on- 
shell, whereas thin are massless. A dot represents 
a squared propagator. 



integrals needed to high orders in the e-expansion 
could be found in the literature. In particular, 
several vertex integrals have been given in |15| . 
However, 7 integrals remained to be evaluated an- 
alytically down to the finite part, Fig. |21 To this 
end, we used series representations in the small 
external momentum regime, which were subse- 
quently improved with the help of conformally 
mapped Pade approximants, and rcsummed em- 
pirically with the PSLQ algorithm. The expan- 
sions were derived from Mellin-Barnes represen- 
tations and analyzed with the help of the MB 
package TB^ . With the usual MS normalization of 



the momentum integrations, i.e. 
per loop, the results are 



h = ^ 
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where S2 = 4/(9V3)Cl2(7r/3), with Cl2(a;) = 
SJ(Li2(e*^)), and the mass has been set to 
unity. The 0{e) parts of these integrals and finite 
parts of the remaining ones have been evaluated 
numerically, either from small momentum expan- 
sions, or from integral representations. There was 
no need to invest into an analytic result. 

During the reduction of the integrals to mas- 
ters, and in particular of a four line propaga- 
tor, an interesting relation has been discovered. 
Fig. 13 Such a relation could certainly not be 
found during the reduction of the three line prop- 
agator graph, since the vacuum graph has the 
same number of lines, but a different mass and 
momentum distribution. In fact, reducing the 
propagator leads to the two masters, the dotted 
and the undotted. This result demonstrates that 
the integration-by-parts relations do not exhaust 
all the rational function coefficient linear relations 
between integrals with a fixed upper number of 
lines in general. 

Note, that the fact that there is a simple re- 
lation between the two 3-line propagator master 
integrals has been known since but there, it 
was not identified as a relation between three in- 
tegrals. 

4. VERTEX CONTRIBUTIONS 

Summing up all of the vertex contributions in 
Feynman gauge, with the mass parameter of di- 
mensional regularization set to A/^, and Mh = 
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Figure 3. Exact relation between 3- line masters [Qgvj 
on-shell. The integral, of which IBPs contain this 
relation is shown in the frame. 



100 GeV, one obtains 

-4.4 X 10-6^(1 - l.ls^ - 3.9s^ - 0.8sV 

+0.2s^ + 1.1 X 10-3 log(s^)s^ + ...) (8) 

Even though just part of the calculation, this for- 
mula is interesting by itself. First of all, the re- 
maining part of the result can be evaluated to 
very high precision for any masses, since it is 
expressed by at most two- loop propagators, and 
therefore, the above series determines the final 
precision. Second, Eq. (|SJ) proves that one should 
resist the temptation to approximate the final 
contribution by the leading term, which is rel- 
atively easily evaluated by setting Mw = Mz in 
all lines. In fact, because ~ 1/4, the sum of 
the first three terms is equal to half of the leading 
one. Finally, the ultrasoft contribution connected 
to the logarithm of sw starts very late in the ex- 
pansion and has a small coefficient, which means 
that one could have obtained a reliable approxi- 
mation without the strategy of regions. 

The convergence of the expansion can also be 
studied by looking at the dependence on M//.The 
first five terms of the expansion are 

- 1.1 X 10-^(1 + O.Usjj + OMs% + 0.02s% + O.Olsf^), 

which shows nice convergence at Mh — 0, corre- 
sponding to = 1. At the other end of the scale 
we obtain an alternating series, which strongly di- 
verges even below the threshold for Higgs boson 



Figure 4. Behavior of the subsequent terms of the 
expansion in of the sum of the vertex graphs. 



decay into a pair of gauge bosons, as shown in 
Fig. El It turns out that this behavior can be of 
advantage. In particular, being alternating this 
series should lend itself nicely to Fade resumma- 
tion techniques, closing the gap between the mass 
difference and large Higgs boson mass expansion. 
A similar phenomenon has been observed during 
the computation of the bosonic corrections to fi 
decay 8^. 

5. SUMMARY 

We have discussed computational techniques 
used in the recent computation of the complete 
bosonic contributions to the effective leptonic 
weak mixing angle. As a by product, we have 
presented new results for on-shell two-loop ver- 
tices and a relation between on-shell propagator 
master integrals. 
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